The off-resonant dielectronic recombination in a collision of an electron with a heavy 

hydrogen-like ion 
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The recombination of an electron with an (initially) hydrogen-like ion is investigated. The effect 
of the electron-electron interaction is treated rigorously to the first order in the parameter 1/Z and 
within the screening-potential approximation to higher orders in 1/Z, with Z being the nuclear 
charge number. The two-electron correction contains the dielectronic-recombination part, which 
contributes to the process not only under the resonance condition for the projectile energy but also 
in the regions far from resonances. The mechanism of the off-resonant dielectronic recombination is 
studied in detail. 
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I. INTRODUCTION 

One of the main processes occurring in collisions of a 
highly charged ion with an electron is radiative recombi- 
nation (RR), in which the electron is captured from the 
continuum into a bound state with emission of a photon. 
In the case when the ion initially possesses one or sev- 
eral electrons, the electron capture can proceed also via 
dielectronic recombination (DR), in which the energy ex- 
cess goes to the excitation of the second electron, which 
then returns to the ground state via a radiative decay. 
DR is a resonant process and is usually studied under 
the condition that the excess energy is very close to the 
excitation energy of the second electron. In this case, 
DR is the dominant recombination channel, whereas RR 
is responsible for a nonresonant background. Outside of 
the resonance region, RR is the dominant process. 

In the zeroth approximation, RR and DR are often 
considered as two independent recombination channels, 
which can be calculated separately and combined addi- 
tively More accurate calculations include the effects 
of quantum interference between DR and RR Q-Q • Gen- 
erally speaking, at the level of precision where effects 
of the electron-electron interaction come into play, RR 
and DR cannot be meaningfully separated. Outside of 
the resonance region, DR is essentially a correction to 
RR due to the electron-electron interaction and induces 
a contribution of the same order of magnitude as other 
two-body effects, e.g., the screening of the nuclear charge 
by core electron(s). 

The accuracy of experimental investigations of the 
RR process with heavy highly charged ions has gradu- 
ally increased during the past years 0-01 1 reaching the 
level on which the electron-electron interaction effect can 
be clearly identified Q. A disagreement with the one- 
electron theory observed in the state-selective study of 
RR into hydrogen-like uranium [sj calls for an accurate 
theoretical description of the electron-electron interac- 
tion effect. 



Most of the previous calculations of the RR process 
into heavy few-electron ions accounted for the electron 
correlation by means of the Dirac-Fock method [1, [l^] , 
disregarding the off-resonant DR mechanism. Evidences 
that the omitted contribution might be significant were 
reported in Refs. [ll|,|l2|, where a part of the off-resonant 
DR (involving photon emission from a core electron) was 
studied. It was claimed that, for many-electron sys- 
tems, this mechanism can significantly influence the RR 
process, yielding an order-of-magnitude enhancement in 
some specific cases. 

In the present investigation we perform an ah initio 
calculation of the electron-electron interaction effect on 
RR into a heavy hydrogen-like ion in the non-resonant 
region of energies of the incoming electron. A particular 
emphasis will be made on the contribution of the off- 
resonant DR, as this effect has not been carefully studied 
before. A similar study of RR into a helium-like uranium 
have been reported previously in Ref. [l3| . 

Relativistic units {h = c = 1) are used in this paper. 



II. GENERAL APPROACH 

We consider RR of an electron with an (initially) 
hydrogen-like ion. The initial state consists of the inci- 
dent electron with the asymptotic momentum p, the en- 
ergy £ = \J + rn^ , and the spin projection fis — ±1/2 
and the bound (core) electron in the state a' with the 
relativistic angular quantum number k^' and the mo- 
mentum projection fia' ■ In the final state, there is the 
two-electron bound state with the angular momentum J 
and the projection M and the outcoming photon with 
the momentum k, and the energy uj — \k\. The wave 
function of the final two-electron state is 
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where a and v stand for the core and the valence electron, 
respectively, and N ^ 1/ \/2 for the equivalent electrons 
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and A'' = 1 otherwise. The core electron state is not 
changed in the process, thus Kq — Kq' . The wave function 
([T]) is not antisymmetrized since we choose to perform 
antisymmetrization exphcitly for the amphtude. 

General formulas are conveniently written in the 
center-of-mass frame, which practically coincides with 
the rest system of the ion. The direction of the z axis 
of the coordinate system is chosen to be the direction of 
the emitted photon. 

In the following, we will assume that the fine-structure 
levels with different J's are not resolved in the experiment 
(as is the case for the experiments conducted so far). 



A. Zeroth order 

To the zeroth order, we neglect the electron-electron 
interaction. The core electron does not participate in 
the process and the transition amplitude is written as 



where t^^)i„ is the amplitude for the recombination with 
the bare nucleus. It reads [141 
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where jw) = denotes the bound state, \p) = |p/is) 

is the Dirac continuum-state wave function with a defi- 
nite asymptotic momentum, and u is the unit polariza- 
tion vector of the emitted photon. After summation over 
the final states and averaging over the initial states, the 
differential cross section of the process is written as 
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where f3 = y^l — vr? je^ . Because of the summation over 
the initial and final states, the cross section does not 
depend on the polarization of the emitted photon, which 
can be fixed arbitrary. The formula ^ differs from the 
corresponding expression for the RR into the bare nucleus 
[l3 | only by a factor of N"^ (=1/2 for the recombination 
into the ground state and 1 otherwise). 

The energy of the emitted photon in Eq. ^ is fixed 
by the energy conservation condition a; = e — or, more 
generally, w = e — m -I- Cio, where Sio is the ionization 
energy of the atom in the final state. 



B. Electron-electron interaction 



correction is induced by a single virtual-photon exchange 
between the electrons, shown diagrammatically in Fig. [T] 
The corresponding correction to the differential cross sec- 
tion can be written as 
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where r^^^^ denotes the first-order correction to the am- 
plitude and is induced by the change of the energy of 
the emitted photon (because of the shift of the energy of 
the final state due to the presence of the second electron). 
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We note that Eq. (O assumes that the perturbative 
regime (|r(^'| <^ I't'-^'I) takes place. 

Since the fine-structure sublevels of the final state are 
not resolved in the experiment, the dependence on J and 
M can be eliminated already in the general formulas. To 
achieve this, we write the correction to the amplitude as 
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where r, 
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does not depend on J and M. Inserting 



this formula and Eq. ^ into Eq. 
summations, wc obtain 



and performing 



dn 



Here, 



da^ 

' dn 



2ja + 1 



4/32e2 



^2Re 



^(0)* _(i) 

'MsMi. Ma Ml. 



r(l) 



MsM. 



(8) 



(9) 



is the amplitude of the recombination with a closed-shell 
atom. So, we obtain that, in the situation when the 
fine-structure levels are not resolved in the experiment, 
formulas for the RR with a hydrogen-like ion differ from 
those for the RR with a helium-like ion only by a prefac- 
tor of N'^/{2ja + l). 

General expressions for the one-photon exchange cor- 
rection to the RR of an electron with a heavy ion were de- 
rived in Ref. [l^. (For the closed-shell ions, such deriva- 
tion was reported also in Ref. ) . The correction to the 
transition amplitude consists of 8 terms corresponding to 
the 8 diagrams on Fig. [U 



For a heavy few-electron ion, the electron-electron in- 
teraction can be effectively accounted for by a pertur- 
bative expansion in the parameter l/Z. The first-order 
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FIG. 1: The one-photon exchange correction to the transition amphtude of the radiative recombination of an electron with a 
hydrogen-like ion. The double line indicates an electron propagating in the field of a nucleus. The wavy line with an arrow 
denotes the emitted photon. The incoming electron is denoted as p, a is the initially bound (core) electron, and v is the 
captured (valence) electron. 



The individual contributions for each diagram are given 
by 
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Here, /(A) is the operator of the electron-electron inter- 
action, 



/(A) = e"a^a.D'"'(A,a;i2) 



(19) 



where D^'^ is the photon propagator. In the Feynman 
gauge, the operator / takes the form 



/(A) ^ ^ l-ai-a2 ^\A\x^, 
47r Xi2 



(20) 



The summations over n in Eqs. (IlT])-([T8]) extend over the 
complete spectrum of the Dirac equation. The second 
term on the right-hand-side of Eq. ([T5|) corresponds to 
the n = V contribution excluded from the summation 
in the first term. The prime on the operator / denotes 
the derivative with respect to the energy argument. The 
state v' is the n = v state with the angular momen- 
tum projection . The small imaginary addition to the 
intermediate-state energies in the energy denominators 
fixes the position of the energy argument of the elec- 
tron propagator G{£) with respect to the branch cuts for 
\£\ > m. 
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We now turn to the physical interpretation of individ- 
ual diagrams in Fig. [TJ The first two graphs represent 
the effect of the screening of the nuclear charge by the 
core electron. The corresponding corrections [t^^'-'^^ and 
T*^^'^^] can be regarded as the first-order perturbations of 
the zeroth-order amplitude ([3]) by the screening potential 
of the core electron, 

where ga and fa are the upper and the lower radial com- 
ponents of the core electron state. The screening effect 
can easily be accounted for to all orders in 1/Z by eval- 
uating the zeroth-order amplitude for an electron in a 
combination of the nuclear and the screening potentials. 
Such treatment is exactly equivalent to the frozen-core 
Dirac-Fock method (as the core in our case contains one 
electron only). 

The contribution of diagram (5) in Fig. [l]can be inter- 
preted to represent the electron correlation on the bound- 
electron wave function (also known as the "relaxation" 
effect). It can be partly included by the standard many- 
body techniques like many-body perturbation theory or 
the multiconfiguration Dirac-Fock method. 

The contribution of diagrams (4) and (8) in Fig. [T] 
contain resonant parts that become prominent when 
the projectile energy approaches the region where e sa 
£v — £a + £n > with e„ being a Dirac bound-state 
energy. When the resonance condition is fulfilled, the 
core electron gets excited into a higher-lying bound state, 
which corresponds to the standard resonant DR mecha- 
nism. In that case, the electron propagator can be re- 
placed by a contribution of the single state responsible 
for the resonance (the so-called "resonance" approxima- 
tion), thus greatly simplifying the problem. In the re- 
gion far from the resonance, however, the core electron 
gets "excited" in all possible virtual states of the energy 
spectrum, so that the usage of the full Dirac propagator 
becomes essential in the description of this process. 

The diagrams (3), (6) and (7) in Fig. [1] correspond to 
other processes with participation of the core electron, in 
which the full energy spectrum of virtual states is probed. 
We will refer to the contribution of all the diagrams (3), 
(4), (6), (7), and (8) as the (off-resonant) DR correction. 
So, in the present work, the term DR is used to refer to 
the recombination with an assistance of the second elec- 
tron, rather than only to the resonant part of this process, 
as is customary. It should be noted that the separation of 
the total two-electron effect in several parts is to a large 
extent artificial (e.g., the DR correction defined in this 
way is not gauge invariant). Its main justification is that 
the screening and correlation parts are easily accounted 
for by standard methods, whereas the DR part is not. 
The sum of all two-electron contributions, however, is 
gauge invariant and derived rigorously within QED. 

So, we represent the total RR cross section as a sum 



of four terms, 

a = + a,er + 411 + 4i : (22) 

where u'-^' is the zeroth-order cross section, asa- is the 
correction induced by the screening potential Vscr in- 
cluded to all orders, Ccorr is the correlation correction 
induced by t^^'^\ and cr^^ is the off-resonant DR contri- 
bution induced by r^^'^), r^i-"), t<^^'^\ t'-^''^\ and t^^'^I 
The screening correction is calculated with the "correct" 
energy of the emitted photon and thus includes the S^i 
correction in Eq. ([S]). We assume that the projectile en- 
ergy is far enough from the resonance condition to ensure 
that the perturbative regime is valid. 

III. NUMERICAL EVALUATION 

The integration over angular variables in the general 
formulas of the previous section can be performed by 
means of the standard Racah algebra, as illustrated in 
Ref. [ll]. The resulting formulas for the zeroth-order 
transition amplitude and for the first-order corrections 
are given in Appendix. Performing our calculations, we 
found several sign mistakes in Ref. [13] . Namely, the 
contributions of Eqs. (|A.8p and (jA.9[) were accounted for 
with the opposite sign in that work. Moreover, the in- 
correct sign was present in the first term of Eq. (|A.6p in 
the case of the capture into the 2pi/2 state. 

The zeroth-order cross section tr*^*^' and the screen- 
ing correction dscr were evaluated according to Eqs. Q 
and (|A.2p . The radial bound and continuum wave func- 
tions were obtained by solving the Dirac equation with 
an extended-nucleus Coulomb potential and the screen- 
ing potential of the core electron, by using the RADIAL 
package by Salvat et al. [16]. 

The calculation of the first-order corrections fJcorr and 
^DR more complicated, due to a larger number of 
radial integrations and the summations over the com- 
plete spectrum of the Dirac equation. In the evaluation 
of the r^^^^^ correction, we ernployed the dual kinetically 
balanced i?-spline basis set [l7| to represent the Dirac 
spectrum. In most of other cases, we used the analytical 
representation of the radial Dirac Coulomb Green func- 
tion in terms of Whittaker functions [l^ . For simplicity, 
we used the point-nucleus Green function, since the ef- 
fect of the finite nuclear size turned out to be negligibly 
small. In the evaluation of the r^^'"'^ and t'^^'^'i correc- 
tions, we used the finite basis set when the energy argu- 
ment of the Green function was smaller than the electron 
rest mass £ < m, and the exact Green function, other- 
wise. The Dirac Coulomb Green function with f > m is 
a complex-valued function and a care must be taken in 
order to choose the appropriate branch of it. The sign 
of the imaginary part of the Green function is fixed by 
the sign of the small imaginary addition in the energy 
denominators of Eqs. (fT^ - P?)l and discussed in detail in 
Ref. QM- 
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A problem emerges in the numerical evaluation of the 
radial integrals when they contain, apart from the Bessel 
function, two continuum-state wave functions. In this 
case, the integrand is a rapidly oscillating function that 
falls off very slowly at large radial distances. In our case, 
such situation arises only in the evaluation of the r^^'^^ 
correction for projectile energies e > m — Sa + £v [The 
problem appears also for the t'^^^^ correction if it is eval- 
uated perturbatively but not if it is evaluated to all or- 
ders.] 

Our scheme of evaluation of radial integrals was as fol- 
lows. First, we introduce the parameter Ri that repre- 
sents the distance at which all bound-state wave func- 
tions become negligibly small. (Typically, i?i = 80 /Z 
a.u.) At the distances r > Ri, all radial integrals with 
bound-state wave functions reach their asymptotic val- 
ues, so that the problem reduces to the evaluation of 
one-dimensional integrals of the form 

poo 

/ drr^ji{ujr)r{r)^^{r), (23) 

JRi 

where ji is a spherical Bessel function, /* is a radial com- 
ponent of the continuum-state Dirac wave function and 
(/)^ is the irregular solution of the Dirac equation (orig- 
inating from the Green function). To evaluate these in- 
tegrals, we introduce a small regulator parameter a > 
and multiply the integrand by exp(— ar). The regular- 
ized integrals are cut off at large distances by a param- 
eter i?2 oc 1/a and evaluated numerically with Gauss- 
Legendre quadratures. The typical value of the regulator 
was a = 10^^. We checked that decreasing the regulator 
by a factor of 10 does not influence our numerical results 
significantly. 

IV. RESULTS AND DISCUSSION 

The calculational results for the total cross section of 
the RR of an electron with an (initially) hydrogen-like 
uranium are presented in Table U for the capture into 
the (Is)^, ls2s, ls2pi/2, and ls2p3/2 states. tT^°^ is the 
zeroth-order cross section. It is calculated with the en- 
ergy of the emitted photon that includes all known one- 
electron corrections to the energy of the final state, i.e., 
Lo = e — m — Sio^Hi where £io,H is the ionization energy 
of the hydrogen-like ion. is the correction due to the 
screening of the nuclear charge by the core electron. It 
was obtained by re-evaluating the zeroth-order cross sec- 
tion with the wave functions calculated in the presence 
of the screening potential. The energy of the emitted 
photon includes all known corrections to the energy of 
the final state jl9|, i.e., lo — e ~ m — ejo, where Sio is 
the ionization energy of the helium-like ion. aiolr is the 
correlation correction induced by r*^^'^). cr^p^ is the off- 
resonant DR correction induced by t^-^'^\ t^-^'^\ t'-'^^^^ 
t'^'^\ and T^-^'^\ For the recombination into the excited 
states, (Tq^ contains a series of the DR resonance peaks 



in the region of projectile energies E = 110—190 MeV/u. 
The behaviour of CTq^ in the vicinity of the peaks is shown 
in Fig. [21 In the case of recombination into the ground 
state, CTp^ does not have any resonances. 

Our calculation shows that the effect of the screening 
of the nuclear charge generally grows for larger projec- 
tile energies and the capture into higher excited states, 
approaching the limit of the complete screening (i.e., the 
case of the capture by a bare nucleus with the Z — 1 
charge). The effect of the off- resonant DR mechanism is 
the strongest for the capture into the ground state and 
for low projectile energies. In this case, the DR contri- 
bution is of the similar size as the contribution of the 
screening effect. We conclude that for the capture into 
the ground state, the electron-electron interaction needs 
to be accounted for rigorously and with inclusion of the 
off-resonant DR mechanism. Results obtained by an ef- 
fective one-electron theory or by standard many-body 
approaches such as the Dirac-Fock method provide only 
an order-of-magnitude estimate of the two-electron effect 
in this case. However, for the recombination into the 
excited states and the projectile energy beyond the DR 
resonance threshold, the DR correction is much smaller 
than the screening contribution and can be neglected for 
most practical purposes. For the projectile energies be- 
low the threshold, the off-resonant DR mechanism can be 
important in the vicinity of the peaks, even at relatively 
large distances from the region of resonance. 

In order to illustrate the dependence of the effects stud- 
ied on the nuclear charge number Z, Table |TT] presents 
the calculational results for the recombination into the 
(Is)'^ and ls2s states of the initially hydrogen-like tin 
{Z — 50). We observe that the relative contribution of 
the screening effect is roughly proportional to 1/Z, as 
could be expected. It is remarkable that the electron 
correlation correction, which plays only a minor role for 
uranium, becomes important for tin in the case of capture 
into the ground state. The relative contribution of the 
off-resonant DR mechanism is slightly larger for tin than 
for uranium, but, in comparison to the screening effect, 
the DR correction becomes somewhat less significant for 
lighter ions. 

In Fig.[3]we present the results for the differential cross 
section for the case of the capture into the ground state of 
uranium, for two values of the projectile energy = 50 
and 300 MeV/u, which are typical for the ESR storage 
ring at GSI. The differential cross section is calculated in 
the laboratory frame, in which the initially free electron 
is at rest. We observe that the screening and the DR con- 
tributions have different dependence on the observation 
angle. For the zero angle, they are of the opposite sign 
and significantly cancel each other, whereas for larger 
angles these two effects amplify each other. 

One of the motivations of the present study was a devi- 
ation from predictions of one-electron theory reported in 
the experimental investigation of RR into a hydrogen-like 
uranium at very small projectile velocities [8]. An effect 
of about 10% was observed in the experiment, whereas 
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FIG. 2: The contribution of dielectronic recombination for the capture into the ls2s, ls2pi/2, and ls2p3/2 states of the 
initially hydrogen-like uranium as a function of the projectile energy E, in units per cent of the zeroth-order cross section a^°\ 
The threshold energies of the resonant dielectronic recombination for the capture into the ls2s, ls2pi/2, and ls2p3/2 states are 
Eo = 178.6, 178.4, and 186.7 MeV/u, respectively. 



a much smaller contribution on the level of 1-2% was 
expected from theory [20j . 

Our ab initio calculation demonstrates that the 
electron-electron interaction affects the RR cross section 
on the level of about 2% for the projectile energies of sev- 
eral MeV/u, which agrees with previous estimates. For 
smaller projectile energies, the cross section is well de- 
scribed by the asymptotic behaviour E a{E) — const, 
and the relative values of all corrections stay constant. 
So, our calculation cannot explain the large two-electron 
effect observed in Ref. i8|. We note, however, that the 
quantities actually measured in this experiment were not 
the cross sections but the recombination rates. A consis- 
tent interpretation of the experimental results requires 
a careful consideration of the recombination rates under 
the experimental conditions. Such a calculation in un- 
derway and will be reported elsewhere. 



V. SUMMARY 

We have performed an investigation of the radia- 
tive recombination of an electron with an (initially) 
hydrogen-like ion. The electron-electron interaction was 
treated rigorously to the first order in the parameter 
\/Z and within the screening-potential approximation to 
the higher orders in 1/Z. The contribution of the off- 
resonant dielectronic recombination was studied in detail. 
It was demonstrated that this mechanism contributes sig- 
nificantly to the total effect of the electron-electron in- 
teraction in the case of recombination into the ground 
state. For the recombination into the excited states, it 
is significant in the vicinity of the resonance peaks but 
becomes small for the projectile energies beyond the res- 
onant dielectronic-recombination threshold. 

The work reported in this paper was supported by 



the Helmholtz Gemeinschaft (Nachwuchsgruppe VH- 
NG-421). 



Appendix: Calculational formulas 

The spherical-wave expansion of the Dirac wave func- 
tion of an incident electron with a fixed asymptotic mo- 
mentum is [13] 

\pHs) = 47r ^ e'^» C^^;^ , Y,*^^ (p) IsKf,) , (A.l) 

where j = \k\ - 1/2, I ^ \k + 1/2| - 1/2, is the 
phase shift, and \eKfi) is the continuum Dirac wave func- 
tion with the relativistic angular quantum number k and 
the angular momentum projection /i, normalized on the 
energy scale. After the integration over the angular vari- 
ables (see Ref. [l^ for details), the result for the zeroth- 
order amplitude is given by 

K. JL 

X V2LTTci^{, i-iy-^ C^Ij-, Pjl{lo,vs) , 

(A.2) 

where the radial integrals Pjl are defined as 

PjL{uJ,ab) = / dxx'^ jLitOx) [gb{x) fa{x) SjL{Kb,-Ka) 

Jo 

~ fb{x) gaix) SjL{-Kb,Ka)] ■ (A. 3) 

The angular coefficients Sjl{i^i, 1^2) are given, e.g., by 
Eqs. (C7)-(C9) of Re f. pTI |. The momentum projections 
/i, Af , and m; in Eq. (IA.2I) are fixed by the selection rules 
of Clebsch-Gordan coefficients. A = ±1 corresponds to 
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FIG. 3: Individual two-electron contributions to the differential cross section of RR into the ground state of the initially 
hydrogen- like uranium, for two values of the projectile energy, iJ = 50 MeV/u (left column) and 300 MeV/u (right column), as 
a function of the observation angle Q. The upper graphs represent the absolute contributions to the cross section in barns/sr 
and the lower graphs, the relative magnitude of the corrections, in units per cent of the zeroth-order cross section da*^"'. The 
dash-doted line (blue on-line) corresponds to the screening part; the dotted line (green on-line), to the correlation correction; 
the dashed line (red on-line), to the DR correction; and the solid line (black on-line), to the total two-electron effect. 



the circular polarization of the emitted photon. (The 
cross section does not depend on the sign of A.) 

The one-photon exchange corrections to the transition 
amplitude r^^'*-* can be expressed in the form similar to 
that for the zeroth-order amplitude, with the radial inte- 
grals Pjh substituted by their generalizations J^^jj^^ ■ The 
results for the functions J^^jf^ are 
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TABLE I: The total cross section of the radiative recombination of an electron into the (Is)^, ls2s, ls2pi/2 and ls2p3/2 states 
of the initially hydrogen-like uranium, for different values of the projectile energy E. a*-"^ is the zeroth-order cross section. (Jscr 
is the screening correction, aioli is the correction due to the electron correlation on the bound electron state, and ctq^ is the 
correction due to the dielectronic recombination. All corrections are given in units per cent of a^'^\ 
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where i?L is the relativistic generahzation of the Slater 
integral (see Appendix C of Ref. (HI). The prime of i?L 
in Eq. (jA.6P denotes the derivative with respect to the 
energy argument, R'j^is, abed) — dj (duj) Rl{<^, 0'^'^'^)\u]=e- 
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TABLE II: The total cross section of the radiative recombination of an electron into the (Is)^ and ls2s states of the initially 
hydrogen-like tin {Z — 50), for different values of the projectile energy E. Notations are the same as in Table |l] 
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